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ABSTRACT 

In order for diffusive shock acceleration (DSA) to accelerate particles to high energies, 
the energetic particles must be able to interact with magnetic turbulence over a broad 
wavelength range. The weakly anisotropic distribution of accelerated particles, i.e., 
cosmic rays (CRs), is believed capable of producing this turbulence in a symbiotic re- 
lationship where the magnetic turbulence required to accelerate the CRs is created by 
the accelerated CRs themselves. In efficient DSA, this wave-particle interaction can be 
strongly nonlinear where CRs modify the plasma flow and the specific mechanisms of 
magnetic field amplification. Resonant interactions have lon g bee n know n to amplify 
magnetic fluctuations on the scale of the CR gyroradius, and lBelll (120041 ) showed that 
the CR current can efficiently amplify magnetic fluctuations with scales smaller than 
the CR gyroradius. Here, we show with a multi-scale, quasi-linear analysis that the 
presence of turbulence with scales shorter than the CR gyroradius enhances the growth 
of modes with scales longer than the gyroradius, at least for particular polarizations. 
We use a mean-field approach to average the equation of motion and the induction 
equation over the ensemble of magnetic field oscillations accounting for the anisotropy 
of relativistic particles on the background plasma. We derive the response of the mag- 
netized CR current on magnetic field fluctuations and show that, in the presence of 
short-scale, Bell-type turbulence, long wavelength modes are amplified. The polariza- 
tion, helicity, and angular dependence of the growth rates are calculated for obliquely 
propagating modes for wavelengths both below and above the CR mean free path. The 
long-wavelength growth rates we estimate for typical supernova remnant parameters 
are sufficiently fast to suggest a fundamental increase in the maximum CR energy a 
given shock can produce. 

Key words: radiation mechanisms: non-thermal — X-rays: ISM — (ISM:) supernova 
remnants — shock waves. 
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1 INTRODUCTION 

Diffusive shock acceleration (DSA) is the most promis- 
ing mechanism for producing superthermal and relativis- 
tic particles in a wide range of astrophysical objects 
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ranging from the Earth bow shock to shocks in galaxy 
clusters jBlandford fc Eichlerl 1 19871 ; Ijones fc Ellison! I Il99ll : 
iMalkov fc Drurvll200ll ~ While this mechanism is believed 
to be efficient and capable of producing cosmic rays (CRs) to 
energies well above 10 15 eV in supernova remnants (SNRs), 
fast and efficient DSA demands that particles interact 
strongly with large amplitude magnetic fluctuations in the 
shock vicinity. The amplitude of the required MHD turbu- 
lence is substantially higher than the ambient MHD tur- 
bulence forcing a bootstrap scenario where the accelerated 
particles generate the turbulence required for their accel- 
eration. Direct support for the modest self-generation of 
MHD turbulence ha s long been seen in heliospheric shocks 
(e.g., iGurnettl [l985h and more recent observations of X- 
ray synchrotron radiation from several young SNRs pro- 
vide indirect evidence for extreme super-adiabatic magnetic 
field amplification (MFA) associated with CR production 
and DSA (see e.g., IVink fc Lamindl2003l : lBamba et al 1 120051 : 
lUchivama et al.ll2007l ; IVinkll2008l ) 

A critical aspect of MFA concerns the dynamic range of 
the self-generated turbulence. Since the maximum CR en- 
ergy a given shock can produce is determined largely by the 
power in the longest wavelength turbulence - that turbu- 
lence which is needed to trap CRs with the largest gyro- 
radii - the production of long-wavelength turbulence must 
be included in any full description of nonlinear DSA. The 
essential features of our calculation for the growth rate of 
turbulence with scales larger than the gyroradii of the gen- 
erating CRs as these CRs interact with turbulence with 
scales shorter than their gyroradii, were prese nted for paral- 
lel propagating shocks in Bvkov et al] l)2009l ). Here, a more 
complete calculation is given which includes oblique shock 
geometry. 

The situation where CRs interact with short- 
wavelength turbulence is just what is expected in a shock 
precursor where the CR current efficiently generates purely 
growing modes with wav elengths m uch shorter than the gy- 
roradii of the CRs (i.e., lBelll2004h . For wavelengths both 
below and above the CR mean free path, we perform a 
multi-scale, linear calculation that takes into account the 
polarization, helicity, and angular dependence of the growth 
rates for obliquely propagating modes. An important ques- 
tion for DSA and CR origin has always centered around the 
maximum particle energy a given shock can produce. For a 
shock of a given size, age, and magnetic field geometry, the 
maximum CR energy depends totally on the power in the 
longest wavelength turbulence. Our estimates suggest that 
the growth rates for long-wavelength modes are fast enough 
to allow a significant increase in CR energy. 

The study of turbulence generation associated with 
CRs and DSA has a long history. Magneto-hydrodynamic 
(MHD) type wave amplification due to the resonant cosmic- 



ray streaming instability was studied in the context of galac- 
tic cosmic-ray origin and 



prop a gation since the 1960s (see 
197ll : IWentzel 1 19741 ; lAchterberd 
"" l2003h . It 



e.g. Kulsrud fc Cesarskyl 

Il98ll ; iBerezinskii et al] Il990l ; IZweibell 12003 ). it was pro- 
posed by iBelll (|l97^ T as a source of magnetic turbulence 
in the test particle DSA scenario, and nonlinear models of 
DSA inclu ding streaming instabilities and MFA wer e inves - 
tigat ed by |Amato fc Blasl (|2006MVladimirov et all (|2006h 
and iReville et all (|2009h . A Monte Carlo model of non- 
linear DSA with MFA from resonant instabilities induced 
by accelerated particles, which also accounted for the ef- 
fects of dissipation of turbulence upstream of a shock and 
th e subsequent precursor plasma heating, was developed 
bv lVladimirov et al] (120081 ) . The Monte Carlo work showed 
that strong feedback effects between the plasma heating due 
to turbulence dissipation and particle injection are impor- 
tant for understanding the nonlinear nature of efficient DSA. 
While the resonant instability is arguably the simplest in- 
stability thus far studied, a full description is still beyond 
any single technique either analytic or computer simulation. 
Mixed techniques are required where analytical recipes to 
account for MFA, dissipation, and other effects are blended 
with simulations. 

In addition to resonant instabilities, a number of non- 
resonant instabilities have been investigated for DSA. A non- 
resonant acoustic instability, where the cosmic-ray pressure 
gradient in the shock precursor amplifes compr e ssiona l dis- 
turbances, wa s inves ti gated by iDorfi fc Drurvl dl985 ) and 
iDrurv fc FaHel ljl986h . iBerezhkol <fl986h and IChalovl jl988h 



generalized this by accounting for a regular magnetic field. 
The obliquely propagating magnetosonic modes in the inho- 
mogeneous precursor s of cosmic- r ay mo dified shocks were 
further examined by IZank et al] (|l990h who also investi- 
gated the role of strong, intermediate, and weak cosmic-ray 
scattering regimes. The effects of the acoustic i nstability on 
the p article distribution were investigated by iKang et al] 
(|l992h using a time-dependent numerical simulations of 
the diffusion-advection transport equations. The non-linear 
evolution of unstable acoustic waves in the precursors of 
cosmic-ray mediat e d sho cks with large Mach numbers was 
shown by IChalovl (|20ld ) to result in the formation of a 
host of small, weak shock waves that could heat the ther- 
mal plasma and, therefore, change the parameters of the 
str ong, large-scale shoc k . Oth er nonlinear work was done 
by IMalkov fc Diamond! |2009f ) where they suggested that 
the development of unstable acoustic waves might result in 
shock-train formation that could speed up the acceleration 
rate and stimulate the inverse cascade of Al fven waves gener- 
ated b y the accelerated particles. Recently. iBeresnvak et all 
(|2009h proposed a model where the stochastic magnetic 
fields in the shock precursor are generated through a short- 
scale dynamo mechanism. In this model, the fluid velocity 
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turbulence driving the dynamo is produced through interac- 
tions of the CR pressure gradient and density perturbations 
in the precursor. 

A large amount of recent work has been stimulated by 
Bell's discovery of a fast, non -resonant instability driven 
by the CR current (|Belll 120041 ) . In this instability, the CR 
current in the shock precursor drives purely growing, in- 
compressible, electromagnetic modes. These modes have 
wavevectors along the magnetic field and have wavelengths 
much shorter than the CR gyroradius radius. Because of 
the short wavelengths, the CR current is essentially unmag- 
netized in this regime and induces a return current in the 
background plasma which is nearly parallel to the locally 
homogeneous magnetic field. This allows the approximation 
that the current generating CRs are only weakly perturbed 
by the magnetic fluctuations they create (we shall discuss 
this instability in some detail in section [3]). 

The focus of much of the analytical and numer- 
ical work on the Bell instability concerns the satura- 
tion level a nd the spectral pro p erties of the instabil- 



ity (see e.g. I Pellet ier et alj l200rj; iMarcowith et al.l l200rj: 



Amato fc B lasi 2009 ;lLuo fc Melrosdl2009l ; IVladimirov et all 



20091 ; IZweibel fc Everett! I2010T ). If the approximation that 



the CR current is only weakly affected by the short- 
scale fluctuations holds, nonlinear MHD-type simulations, 
where the CR cu rrent is fixed as an external parameter , 
can be used (e.g . [Belli l200l; IZirakashvili fc Ptuskinl l200d: 



IZirakashvili et all |200S| ; iReville etldl l200ct ). These simula- 
tions typically show a rather high saturation level of the 
self-generated, short-scale turbulence and also show nonlin- 
ear spectral energy transfer (i.e., the cascading of turbu- 
lence energy) into both larger and smaller scales. However, 
since long-wavelength magnetic fluctuations on scales com- 
parable to or larger than the CR gyroradius cause a non- 
negligible response of the CR current, the MHD simulations 
are limited to short-scale fluctuations only. In a real system, 
long-wavelength fluctuations will occur and they will induce 
current perturbations perpendicular to t he local mean mag- 
netic field even in the li near regime (see iBvkov fc Toptyginl 
l2005l ; lBvkov et al.ll2009T ) . This important physical effect will 
be difficult for MHD simulations to capture. 

Particle-in-cell (PIC) and hybrid plasma simulations^ 
are widely used with different CR and background plasma 
parameters. They can explore both linear and nonlinear 
regimes of the non-resonant mode, the sa turation level, 
and s pectral energy transfer mecha ni sms (iNiemiec et al.l 
20081; iRkmehne fc Spitkovskvl 120091; lOhira et all I 2OO9F 



Stroman et"ai] 120091 ; IRiguelme fc Spitkovskvl I201C 



1 In this context, hybrid means protons are treated as particles 
and electrons are modeled as a charge-neutralizing background 
fluid. 



IGargate et all l2010h . In the linear regime, the simu- 
lat ed gr o wth r ates confirm the basic theoretical predictions 
of IBelll (|2004f ). Amplification factors above 10 for the 
magnetic field were achieved in some simulations. Plasma 
acceleration along the drift motion of CRs saturates the 
instability at the magnetic field level when the Alfven 
velocity in the amplified field approaches the drift ve- 
locity of CRs. The instability can also saturate earlier 
if CRs get magnetized by the amplified field. If the CR 
current is transverse to the initi al magnetic field then both 
hydrodynamic (e.g., Bell 20051 ) and PIC simulations by 
IRiguelme fc Spitkovskvl (|2010l ) have shown the presence of 
a fast, current driven instability. 

Magnetic field amplification in the long-wavelength 
regime, where the gyroradii of the "magnetized" CRs are 
smaller than the wav elength of t h e ma gnetic perturba- 
tions, was discussed bv lBvkov et al.l (|2009l ). This work used 
a quasi-linear kinetic equation approach for modes paral- 
lel to the in i tial m agnetic field. It was demonstrated by 
IBvkov et all (|2009l ) that the long-wavelength instability, 
produced by the perturbed, perpendicular CR current, was 
strongly affected by the short-scale turbulence produced by 
the Bell instability. This coupling of long and short-scale 
instabilities, which is certain to be nonlinear if DSA is ef- 
ficient, places significant demands on plasma particle and 
fluid code simulations since they must cover a wide dynamic 
range to fully describe the problem from injection to maxi- 
mum CR energy. These demands are particularly severe for 
nonrelativistic shocks w ith parameters typical of SNRs (see 
IVladimirov et al.ll2008l ). 

In the present paper we study the properties of long- 
wavelength, oblique modes in DSA for different assump- 
tions on the CR scattering rate in background, short-scale 
turbulence. The long-wavelength growth rates we calcu- 
late depend on the CR current interacting with the short- 
wavelength turbulence generated by Bell's instability, which 
has the fastest growth rate in the short wavelength regime 
and produces the mode polarizations we assume. While 
our calculation is linear in AB/B, any application in a 
realistic DSA scenario would assume efficient acceleration 
with a non-negligible fraction of shock ram pressure being 
transformed to CRs. Our estimates given in Section [6] for 
young SNRs, suggest that, for some parameters at least, 
the production of long-wavelength turbulence in the for- 
ward shock precursors of young SNRs will be important. 
Since the shock structure depends critically on the efficiency 
of turbulence amplification, cascading, and dissipation, our 
work must be viewed as a step towards the development 
of a nonlinear model of DSA where all of the following are 
treated self-consistently: (1) efficient particle injection and 
acceleration occur; (2) particles of different energies par- 
ticipate differently in the instability generation and mag- 
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netic field amplification; (3) turbulence cascading and dissi- 
pation are accounted for; and (4) the nonlinear feedback 
of particles and fields on the bulk flow is included self- 
consistently. To our knowledge, the only techniques that 
are currently capable of including all of these coupled ef- 
fects in a calculation suitable for modeling DSA in SNRs, 
even in parameterized form , are the semi-analytic methods 
of Bl asi and co-workers (see iBlasi et al.ll2007l ; ICaprioli et alj 
120081 . and references therein) , and t he Monte Carlo methods 
of Vladimirov and co-workers (see IVladimirov et all |2009| . 
and references therein). Before providing the details of our 
derivation we outline the steps in the next Section. 



2 OUTLINE OF THE MODEL 

Our goal is to obtain the growth rates for the cosmic-ray 
current driven, long- wave instability in a shock precursor 
containing short-scale fluctuations. The terms "long" and 
"short" are relative to the minimum CR gyroradiuqj and 
the first step in our derivation is to obtain the equations of 
the plasma and magnetic field dynamics averaged over the 
short-scale fluctuations. 

In Section[3]we describe the properties of the short- scale 
CR-current driven modes produced by Bell's instability (|Belll 
2004, 2005). Wc discuss the polarization properties of the 
modes, the helicity they induce in the background plasma, 
and emphasize how the mode polarization is important for 
the long-wavelength dynamics. 

The equation of motion and the magnetic induction 
equation, both for the background plasma and both aver- 
aged over the short-scale modes that are needed to study 
the long-wavelength growth rates, are discussed in Section 
[4] The detailed derivations of these averaged equations are 
given in the Appendix. In Section [C] of the Appendix we 
obtain the averaged magnetic induction equation. In Sec- 
tion [D] we average the momentum equation. We general- 
ized th«_jneaii_^«ld_meth£d developed in dynamo theory 
fsee iBlackman fc Field! |2002| ; Brandenb urg fc Subramanianl 
120051 ; lBrandenburgH2009T ) to average the equations account- 
ing for the specific effects of the CRs. The distinctive feature 
of the derivation is that the CR current, which is responsi- 
ble for the instability, is fully accounted for. The CR current 
produces specific polarizations of the modes and this influ- 
ences the ratio of their kinetic and magnetic energy densities. 
This, in turn, modifies the mode correlation functions that 

2 Due to the strong energy dependence of the CR diffusion coef- 
ficient, the CR distribution in a shock precursor typically shows a 
fairly sharp cutoff at low energies. This is the minimum energy we 
refer to and this cutoff energy increases with distance upstream 
from the viscous subshock. 



determine the turbulent transport coefficients in the mean 
field approach. The correlation functions are presented in 
Section |B] of the Appendix. 

To get the linear dispersion relation for the long- 
wavelength fluctuations, one needs to know the response 
of the CR current to the short-scale fluctuating magnetic 
field, that is, to the Bell turbulence. This is an important 
point. For solely the Bell instability, the response of the CR 
current to the short-scale fluctuations can be ignored. In the 
long-wavelength regime, however, the CR current changes in 
response to the magnetic field perturbations imposed on the 
plasma system and these changes need to be accounted for. 
Here, we calculate the CR response using a kinetic equation 
for the CR distribution function with a collision operator 
that describes the CR scattering by magnetic fluctuations. 
The details of the CR current response derivation are pre- 
sented in Section [X] of the Appendix. 

Finally, the growth rates we obtain for long-wavelength 
fluctuations in plasma systems with a CR current are pre- 
sented in Section [5] Two distinct long-wavelength regimes 
are discussed. The first regime is where the growing modes, 
at some position in the shock precursor, have wavelengths 
between the gyroradius and the mean free path of the lowest 
energy CRs present at that precursor position. We call this 
regime "intermediate" and these modes are complementary 
to those seen in mean-field dynamo theory except now the 
modes are modified by the presence of the CR current. These 
long- wavelength modes are produced by the anisotropic CR 
distribution in the presence of strong, small-scale fluctua- 
tions. The second regime is where the growing modes have 
wavelengths longer than the CR mean free path. We refer to 
this regime as "hydrodynamical." We note that in the case 
of very strong CR scattering, that is when the Bohm limit is 
obtained and the CR mean free path equals the gyroradius, 
the intermediate regime has no dynamic range and only the 
hydrodyamical regime produces the long-wavelength insta- 
bility. 



3 SHORT-SCALE DYNAMICS AND BELL'S 
INSTABILITY 

In MHD-type flows with cosmic rays imbedded in a back- 
ground plasma, the momentum equation of the background 
plasma, including the Lorentz force, is given by 

p(j^- + (uV)uj = - VP + i (j x B) + e(m - n e )E + vAu, 

(1) 

where u and j are the bulk velocity and the electric current of 
the background plasma, respectively. The viscosity v is due 
to Coulomb collisions (or due to plasma oscillations on scales 
much less than the CR gyroradii) and v is typically small 
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for the effects discussed here. We assume quasi-neutrality 
for the whole system consisting of background plasma ions 
of number density m, electrons of number density n e , and 
cosmic rays of number density n cr . For simplicity we con- 
sider cosmic-ray protons only such that m + n cr = n e . Both 
the background electric current j and the electric current of 
accelerated particles J cr are the sources of magnetic fields 
in Maxwell's equations, where the displacement current can 
be omitted for the slow MHD-type processes under consid- 
eration. Furthermore, we assume ideal plasma conductivity 
in Ohm's law: 



field and 



rotB 



4-7T ,. 



J"'), E 



--(u x B) 



Then the induction equation is given by 



<9B 
~dt 



V x (u x B) + u m AB, 



(2) 



(•3) 



where v m is the magnetic diffusivity due to the Coulomb 
collisions or MHD plasma oscillations. 

To study instabilities in flows with initially quasi- 
homogeneous magnetic fields, the global current neutraliza- 
tion condition should be fulfilled in Eq. ([2]). Substituting 
Eq. ((2j) in Eq. (fT|). we obta in the momentum equation in the 
following form (Bell 2004) 



p(| + (uV)uj =-VP+i(VxB)xB- 
-i(J cr - en cr u) x B + i/Au . 



(4) 



In this equation, the electric current of accelerated par- 
ticles J cr is an external current in the momentum equa- 
tion of the background plasma. The current is governed by 
sources of energetic particles and by local electromagnetic 
fields. The current J cr initiates a compensating return cur- 
rent in the background plasma. Bell ( 2004, 2005) discovered 
that the system is unstable against linear perturbations that 
are oc exp{^t + ikr) and the return current drives nearly 
purely growing modes. Here, 7 is the linear growth rate. 
The wavenumbers of growing modes must satisfy the condi- 
tion kr g o > 1, where r 9 o = cpo/(eB^) is the gyroradius of 
an accelerated particle of momentum po- All particles with 
p > po contribute to mode growth. 

In a cold plasma with sound speed ao, well belo w the 
Alfven velocity v a , the linear growth rate obtained by iBelll 
|2005h depends only on the wavevector projection k z on the 
local mean magnetic field, i.e., 



7 ~ 7max^/^i 



(5) 



where 



7max = Va \/fcl|fcj — fc 2 (6) 

is the growth rate for the modes propagating along the mean 



fci = 



4-7T j c 



B 



(7) 



Here, the bar means the CR current and magnetic field are 
averaged over fluctuations with scales bel ow r„ Q. 

According to the linear analysis of IBelll (|2005h . the 
wavenumber of a growing mode must satisfy the condition 
r g0 < k < ki. Therefore, the instability growth rate is 
higher than the Alfven frequency v a k. Note that this con- 
dition for mode growth r g ofci > 1, together with Eq. (0, 
implies that the anisotropy of the relativistic particle dis- 
tribution, 5 cr , exceeds the ratio of the mean magnetic field 
energy density to the energetic particle energy density E CI . 
That is, Scr > B 2 /(4:ivE CT ), where the CR current is given 
by j CI « 8 CT en cr c. 

The polarization of the growing mode, given by 



bx=l h bv > 



(8) 



is important and we use the fact that the kinetic energy 
density in the growing mode dominates over the magnetic 
energy density to get simplified mean field equations. The 
linear relation between the amplitude of the growing mag- 
netic field b(k) and the velocity v(k) of the perturbations 
is given by 



This yields 



7 b(k) = iB k z v(k). 



'^' 2 -4^ |b(k)|2 ' 



(9) 



(10) 



provided that the kinetic energy density in the growing mode 
dominates over the magnetic energy density because ki > 
k z . This is in contrast to Alfven modes where the energy 
densities are equal. In Section [B] [namely, Eq. (|B2[l ]. the 
mode polarizations are used to express the pair correlation 
functions for the fluctuating fields. 

We note here that the cosmic-ray current has only a 
weak response to the short-scale fluctuations, while its re- 
sponse to fluctuations with scales longer than r g o is much 
larger. Therefore, the CR current variations can be neglected 
when averaging the dynamic equations over fluctuations pro- 
duced by the Bell instability and we perform this averaging 
in the next Section. For the long-wavelength fluctuations, 
the current variations must be considered and this averag- 
ing is done in Section A of the Appendix. 
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4 THE AVERAGED EQUATIONS OF 
LARGE-SCALE DYNAMICS 

Here we average the momentum equation of the background 
plasma and the induction equation over the short-scale fluc- 
tuations produced by Bell's instability to obtain the mean 
field dynamics equation. Since the linear growth rate is fast 
for the Bell mode with the wavevector along the local mean 
magnetic field, the bulk velocity of the background plasma 
u(r, t) and B(r, t) both experience rapid, incoherent fluc- 
tuations on scales that are small compared to r 9 o, the gy- 
roradius of a relativistic particle. The unstable modes with 
wavevectors along the mean magnetic field are nearly in- 
compressible making the averaging of the mass continuity 
equation straightforward. The time and spatial scales sepa- 
ration relations are defined by J cr = j cr +j cl , B = B + b and 
u = V + v, where j cr , B, and V are the averaged electric 
current of accelerated particles, the averaged magnetic field, 
and the background plasma bulk velocity, while v and b 
are the short-scale bulk velocity and magnetic field. We can 
now average over the ensemble of short-scale fluctuations, 
placing these averages in angular brackets, and obtain the 
averaged momentum equation Eq. (Q: 

^ + (VV)V = - <(vV)v> + ^((Vxb)xb)- 

-IvP- — (Q^-en CT V) x B) + 
P CP 



+ ^((VxB)xB), 
and the averaged equation of magnetic induction 



&B 

dt 



cV x 8 + V x (V x B) 



■AB. 



(11) 



(12) 



Here £ — (v x b) is the average turbulent electromotive 
force. The coordinate system is defined relative to the un- 
perturbed magnetic field e z = B /-Bo- 



5 THE LONG- WAVELENGTH COSMIC-RAY 
CURRENT DRIVEN MODES 

In the presence of a cosmic-ray current, Bell's instability 
results in the fast growth of short-scale modes with wave- 
lengths shorter than the gyroradius of the cosmic-ray par- 
ticles. However, as we showed above when we obtained the 
mean field dynamic equations averaged over the ensemble 
of short-scale motions, the strong short-scale turbulence in- 
fluences the plasma dynamics on scales larger than the CR 
gyroradii producing the turbulence. Equations (|f II) and (|12|) 
are designed to be applied to the dynamics of modes with 
scales larger than r g o, i.e., CR particles are magnetized on 
these scales. 



In diffusive shock acceleration with strong scattering, 
the particle mean free path, A, is often taken to be A = 7]r g o, 
where r\ I and rj = 1 is the Bohm limit. If r\ > I there 
are two different regimes for the large-scale dynamics. The 
first regime, discussed next, corresponds to rj^ 1 < kr g o < 1, 
where the CR particles have gyroradii small compared to 
the turbulence scales and can be considered magnetized. The 
second regime is for kr g o < g^ 1 , where the particle gyroradii 
are large compared to the turbulence wavelengths. These 
modes are driven by the transverse current of anisotropic, 
magnetized cosmic rays and are discussed in Section [5[2] 



5.1 Long- wavelength current driven modes in the 
intermediate regime (g^ 1 < kr g o < I) 

The dynamic equations averaged over the short-scale fluc- 
tuations, that we obtained in Sections [C] and [D] can now 
be used to study the long-wavelength modes. We start with 
the linearized equation of motion Eq. (|D5[) an d the induc- 
tion equation for the mean field Eq. (|C7[) . Denoting the small 
departures of physical values from their unperturbed mag- 
nitudes by 8, and performing the analysis in the rest frame 
of the unperturbed upstream plasma so that V = <5V, we 
obtain: 

9SV 1 „ 1 ,. cr 

-— = — VP (jo e* x <5b - 

ot p cp 

-—(Sj cr x Bo e z ) + J_((V X 5b) x See*) + 

Cp VKp 

+ ~cp ( ~ 9 ' Sbx ^ Gy ~ ~ 9 >Sby ^ ' ^ 



35\) 2 — cr 

— — = V x (SV x B e z ) + —a t V x j e 2 + 
ot pc 



+ 2pc" tV X ^^ ex + 5 ^y e y) + 
3 

+ - — atV x (g'Sb x e x + g'Sb y e y ) , 
Zpc 



(14) 



where g' — jjf/Bo = en CI v s / Bo, v a is the shock velocity, 
and at = (b 2 )T cor . In these linearized equations we have ne- 
glected terms ~ en cr SV compared to g'Sh which is justified 
if Vph/v 3 <C 1 (where v p h = u/k is the phase velocity of the 
modes defined in Section IA1[) . The correlation time r cor is 
the relaxation time of triple correlations (see Sections [Cl and 
[Dl in Appendix, and namely Eq. IC4[) and is approximately 
the turnover time of the Bell turbulence. The parameter r cor 
is an important parameter in the m ean field approach (e.g., 
iBrandenburg fc Subramanianll2005h . 

From Eqs. (|13I14[) . and using Eq. (|A10[) . one obtains 
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the dispersion relation in the form: 



2 -r- -11 Qt 

47rp 



±kk v't ( 1 + ) [A(a!o) - 1] = 



1 „/ n 3 



(15) 



where fco = 4irg'/c, v a = Bo/y/Anp, the function A(xo) is 
defined by Eq. (|A9[) m Section IAU and Kt is a turbulent 
transport coefficient in the averaged equation of motion and 
is defined in Eq. (|D5I) . The ± signs give the modes with the 
opposite circular polarizations in the dispersion relation. 
To estimate the coefficients in Eq. (|15[) we assume 



(b 2 ) « Bq and introduce a dimensionless parameter for the 



amplitude of the Bell turbulence Nb = 



2 )/B . The tur- 



bulence correlation time r cor , is estimated as the short-scale 
mode vortex turnover time, while the amplitude of the tur- 
bulent velocity is (v 2 ). Then, the turbulent mixing length 
is defined as r CO r\/ (v 2 ) w 



! ) w r cor ^(b 2 )/(47rp) = 2^/fc 
where £ is a dimensionless parameter characterizing the tur- 
bulence mixing length. In our numerical estimations below, 
we assume £ ~ 5 for the short-scale turbulence produced by 

Bell's instability. 

According to lBelll |2004l ). the maximum growth rate of 
the instability is at the wavelength A = 47r/fco. However, due 
to the nonlinear dynamics of the Bell instability, the maxi- 
mum of the turbulent energy density may be at somewhat 
larger wavelengths if the initial spectrum of the fluctuations 
subject to growth is a declining function of the wavenum- 
ber k (c.f.. IZirakashvilFet ah 2008; iRiauelme fc Spitkovskvl 
120091 ). We specify the scalings of the turbulent kinetic coef- 
ficients with the dimensionless parameter £ as fcoQt/(47rp) ~ 
2-k^NbVo,, Kt/Bo ~ tvNb, and assume that the minimum 
momentum of accelerated protons is po — mc. With these 
values, the solution to the dispersion equation Eq. (|15[) is 



4c 



with 



and 



=fikk -^- 

47rp 



-k 2 Va ± kkoVa 



1 a, ^ 3 

- 2 A(xo) + - 



1 + 



Kj 

Bo 



[A{x ) - 1] . 



(16) 



(17) 



(18) 



Below, we present simple estimations for frequencies, 
growing-rates and polarizations for the fastest growing 
modes. As ha s long been known, the standard resonan t in- 
stability (e.g- lBlandford fc Eichle3ll987l : lKulsrudll2005h op- 
erates in the intermediate magnetized regime ?7 _1 < kr g o < 
1 in the lack of strong, short-scale turbulence. This reso- 
nant effect in the dispersion equation (|15|l is dominated by 
the imaginary part of the current response function A(xq) 



defined by Eqs. (|A9l IA121) . The response function was cal- 
culated for a power-law momentum distribution of CRs of 
index a as defined in Eq. IA2I Modes with different circu- 
lar polarizations that are distinguished by the sign =p in 
Eqs. (|A9IIA12|) have the same growth rate and with no short- 
scale turbulence (i.e., at — and nt = 0), the frequencies of 
the two circularly polarized modes are determined by 



±(1 + 



3tt 



1 



kor g o(kr g0 



(a-2)/2 



(19) 

This growth rate is plotted as the dashed curve in Fig. [Tjfor 
kor g o = 100 and a — 4. 

Now let's consider the effect of strong, short-scale Bell 
turbulence on the resonant instability when at and Kt are 
nonzero. The main contribution to the dispersion equa- 
tion (|16|l at £ ~ 5 is due to the coefficient d in Eq. (|17p . 
Since the response function A(xo) ~ 1 for kr g o < 1, the 
growth rate can be approximated as 



= i4TT\flN B Vak 



(20) 



For kr g o < 1, only the mode with polarization Sh = 8b{e x + 
ie y ) is growing while for the case kr g o > 1, Bell's instability 
amplifies the other mode with polarization Sh = Sb(e x — 
ie y ). The ratio of the kinetic energy density to the magnetic 
energy density in the growing mode can be estimated from 



|*V(k)|' 



47T/90 



3fc r 9 o(l + 7r7V B ; 
Mtv^Nb 



\6b(k)\ 



(21) 



It should be noted that the helicity of the unstable, 
long-wavelength mode studied above is opposite that of the 
short-scale Bell mode. This provides, in principle at least, 
the possibility of balancing the global helicity of the system 
by combining short and long-wavelength modes. Care must 
be taken however, since recent numerical simulations show a 
high saturation amplitude of the Bell mode making a nonlin- 
ear analysis necessary to address the helicity balance issue. 
The estimations given above are valid in the intermediate 
regime and provide simple analytical approximations to the 
growth rates shown in Fig.[T]for kr g o ^ 1. To turn to the hy- 
drodyamical regime kr g o < rj^ 1 (considered in Section [5. 2[) , 
one should just change (1 + m/Bo) to K, t /Bo in Eq. (fT8)) . so 
that 



1 2 2 , 
-k v n ± 



kk v a K t [A(x ) - 1] /B 



(22) 



and then substitute this in Eq. (|16[) . 

The coefficient given by Eq. (|22[) (or more exactly the 
imaginary part of Eq. IA11|) dominates the dispersion equa- 
tion Eq. (|16p if xo <S 1 and 77 is finite. Because of the =F 
sign in the imaginary part of Eq. (|A11[) . both circular po- 
larizations will be growing with the same growth rate given 
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by 



7 ! 



tvN e 



kko 




n 



(23) 



The transition from the regime described by Eq. (|20l) to that 
described by Eq. (|23[) takes place at 



Xq ~ -- 



?7 32ir(y/NB 



(24) 



5.2 Long- wavelength current driven modes in the 
hydrodynamical regime (kr g o < rj^ 1 ) 

The ponderomotive force ((j cr — en cr v) x b)/(cp) in the 
mean-field momentum equation of the background plasma 
Eq. JD5J, is due to the momentum exchange between the 
background plasma and cosmic rays. Contrary to the short- 
wavelength regime, the cosmic-ray current response on the 
magnetic fluctuations is essential in the long-wavelength 
regime, kr g o < 1, and results in a non-negligible pondero- 
motive force. If the perturbation wavelength is longer than 
the CR mean free path, A, or in other words if kr g o < J?" 1 , 
then the hydrodynamic approximation can also be applied 
to the cosmic-ray dynamics. Then, the momentum density 
P' r ' and the momentum density flux tensor II ^ of the CR- 
fluid, defined as 



n S = j ^P/3/( r 'P' t ) d3 P = PrU^U^ + P 



cr<J a/ 9, 



(25) 



can be approximately derived in closed form. Here /(r, p, t) 
is the CR distribution function. Then, the CR momentum 
equation, derived as a moment of the kinetic equation, takes 
the form 



at 



+ v fl n 



(r) 



-[(J cr X B) + en cr E] Q . 



(26) 



Equation (|26[1 , averaged over the fluctuations with scales be- 
low the CR mean free path and taking into account Eq. ([2]), 
has the form 



OP 



0) 



Of 



-[O cr - en cr V) xB + ((f- en cr v) xb)]„. (27) 



The CR distribution in Eq. (|25[) is nearly isotropic for 
scales larger than the CR mean free path kr g o < Then, 

(r) 

the isotropic cosmic-ray pressure dominates in 111^ and us- 
ing Eq. (|27|) one may write a simplified Eq. (|D5f) in the form 



= _i V (SP + 8P cr ) + J-((V X Sb) x B e z ) + 
ot p inp 



cp 



■ g'Sb x ] e y - [8fy - g'Sb y ] e x ) 



(28) 



In the rest frame of the unperturbed shock precursor, the 
velocity has only the perturbed component V = <5V. The 
term ~ en CI SW is small compared with g'Sh if v p h/v s <C 1, 
and it is omitted in the dispersion relation. 

To study the dispersion properties of the CR current 
driven modes in the long-wavelength regime, we use the 
mean-field dynamic equations averaged over short-scale fluc- 
tuations, Eqs. ()C7[1 . and (|28fl . and the mass continuity equa- 
tion: 



dp 
dt 



+ V(pV) = 



(29) 



The equation of state for the background plasma is assumed 
to be adiabatic, i.e., 



dP 
~dt 



+ (VV) P + 7 PVV = 0, 



(30) 



where 7 is the adiabatic index. Using Eq. (|29p and Eq. pop , 
one obtains 



V5P = atVSp, 



(31) 



where ao = \J 7P0 Jpo is the sound speed of the background 
plasma. 

To get the dispersion equation one must substitute the 
expressions <5b = B — Bo, <5V = V, <5j = j cr — jy into 
Eqs. (EH and (gSJ). Then, using Eqs. (jAT9|) and J3TJ, by 
neglecting the terms that are small at x < xo (where xo is 
defined by Eq. I24[) . we finally get the linear dispersion rela- 
tion for the perturbations with wavelengths larger than the 
mean free path of the cosmic-ray particle A, i.e., 

(u) 2 — k 2 v 2 — ik\\BoWo){uj 4 — (a 2 , + vl)k 2 oj 2 + 

+k\k 2 v 2 a 2 l — ifc||B [(^ 2 — a 2 k 2 )W^ + k±k^alW4]} — 

-fcfBoWiKoj 2 - k\a 2 )W2 - k x k^alW 6 ] = 0, (32) 

where we assumed k x = 0, k y = k± , and 

Wo = g' — -Di(k h k x ), 
cpo V 



Wl=g '^L [D^kx)-!] 
cp 



W 2 



g'— \D 1 (k l i,k A 
cp 



1 



,_Bp_ 1 ^ ^ 

cpo V 



W 3 = g' — -[Dx(k h k±) + Da(*||,fcx)] + 
cpo V 

+g'^l D2 (k h k ± ), 

cpo ?7 

WA = g'—D 3 {k h kx), 
cpo 

W 5 =g'—^D 3 (k h k±). 

cpo n 



(33) 



The angular dependence of the coefficients Wi for i= 0, 
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4 is determined by 
£>i(ft||,fci_) 

D 3 (k h k±) 



k-'i 



fcjj + 
k i k 



p-, _D 2 (fc||,fc x ) 



t.2 , fc 2 ' 



S.2 , fc 2 



(34) 



All of the coefficients Wi given above are proportional to 
g' and thus they vanish in the absence of the unperturbed 
cosmic-ray current, i.e., if j" — 0. The corresponding terms 
are therefore responsible for the current driven modes. The 
coefficients have different asymptotic behaviors at r\ 3> f. 
The coefficients Wo, W3, and W& scale oc r)^ 1 , while all of 
the others are oc r)~ 2 . It is important to note, however, that 
contrary to W3, W4 oc k\\ and Wo oc fejf. 

The growth rates depend on the mode propagation an- 
gle and we first consider parallel propagating modes. As dis- 
cussed above, for modes propagating parallel to the initial 
magnetic field, the growth rate is 



70 (k) 




(35) 



and these modes have the fastest growth rates for the Bohm 
diffusion regime with r\ ~ I. In Fig. [T] we illustrate the effect 
of short-scale turbulence on the long-wavelength instability 
for a particular set of parameters. The dashed and dotted 
curves show the result without short-scale turbulence. The 
two curves show two growing short-scale modes, two other 
modes in Bell's dispersion equation are not growing and 
aren't shown. The departure of the solid and dot-dashed 
curves from the linear dependence (evident flattening) at 
small k is due to the transition to the hydrodyamical regime 
when kr g o < rj^ 1 . It is clearly seen in Fig.[T]that the resonant 
growth rates calculated for the model with short-scale tur- 
bulence are well above the standard resonant growth rates 
shown by the dashed line at kr g o < 1. 

If 77 > 1, scattering is less efficient than the Bohm 
limit and the maximum growth rate occurs for obliquely 
propagating modes. The non-parallel propagating modes are 
driven by both the cosmic-ray current and the cosmic-ray 
pressure gradient (i.e., the diffusive part of the cosmic ray 
current). An analysis of the relative contributions of the cor- 
responding terms in W3 defined in Eq. ((33)) using Eq. (|AI9|) . 
shows that, at the maximum growth rate, the relative con- 
tribution of the cosmic-ray pressure gradient to that of the 
cosmic-ray current in the shock precursor is oc ttNb ■ There- 
fore, the cosmic-ray current contribution is larger if the 
short-scale Bell turbulence is strong. In Fig. [2l the angular 
dependence of the hydrodyamical, long-wavelength modes 
is illustrated for a finite temperature plasma with the pa- 
rameter /3 = ao/^a = 1-0. One can see that, in contrast to 
the short-scale Bell instability, the modes propagating along 




1(T 3 1(T 2 10 4 1.0 10 



kr 



so 



Figure 1. The figure shows growth rates of the parallel propa- 
gating modes as a function of the wavenumber to illustrate the 
effect of short-scale turbulence on the long-wavelength instability. 
Equation 11161 was solved numerically to generate the curves. The 
model parameters are fcor 9 o = 100 and a = 4.0. The solid and 
dot-dashed curves are simulated for two modes in the model with 
short-scale turbulence of £ = 5 and rj = 10 to demonstrate the 
behavior of the modes in the intermediate regime. For compari- 
son, the dashed and dotted curves are calculated for the model 
with out the sho rt-scale turbulence, i.e., with Ng = and r\ 00 
fc.f., lBellll2004r) . 



the unperturbed magnetic field are growing while the modes 
propagating in the opposite direction are damped. It is also 
clear in Fig. [2] that the fastest growing modes are propa- 
gating nearly perpendicular to the unperturbed field Bo for 
r\ — 10, a typical result for 77 2> I. 

To get the propagation angle <9 ma x of the mode of max- 
imum growth for r\ > 1 , one needs to find the maximum of 
the expression k^Ws(k^/k = cos#), i.e., 



cos 9 n 



1/V, 



(36) 



and the maximum growth rate at 77 3> 1 is determined by 



(37) 



The results were obtained assuming v p h/v 3 « 1. In our 
case 7(fc) ~ w(fc), and therefore from Eq. 1371 one may get 
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Figure 2. The two panels show the angular dependence of 
the growth rates of the long-wavelength unstable modes in the 
hydrodyamical regime where kr g o < r]^ 1 . In the left panel 
r] = 1, while r\ = 10 in the right panel. The model param- 
eters are kor g o = 100, a, = 4.0, and the plasma parameter 
= Oq/v^ = 1.0. The solid and dot-dashed curves show the two 
unstable modes. The normalizing parameter 70 is determined by 
Eq. g5}. 



« 1 



the validity condition in the form 

(where M a is the Alfven Mach number of the shock) . 

Consider the polarization of a mode propagating along 
the direction giving the maximum growth rate for rj ^> 1 
(i.e., k± 3> fcii, see Fig. [5J). The amplitude of the magnetic 
field along the initial magnetic field, 8b z , in this mode, ex- 



ceeds the transverse magnetic field perturbations 8b x>y . The 
maximum velocity component SV y is in the plane deter- 
mined by the wave vector and the initial magnetic field 



5b- 



5b v 



SV = - 

V V a k± y/4npo v a k\\ y/Airpo ' 



(38) 



and the energy density of the mode is dominated by the 
kinetic energy 



\SV y (k)\ 2 



47rpo 4 



\8bz(k)\ , 



(39) 



since fco S> k, and we assumed |fej_| ~ k. 

The modes are compressive and from the continuity 



equation one can write 

8_p 

Po 



k ± 5V y + k {] SV z 



k ± 8V v 



(40) 



Sp 




Sb 


Po 




Bo" 



Then, using the induction equation given above, the equa- 
tion V<5b = 0, k± > fey, and 8V y 3> 8V Z , one can estimate 

The complex frequency oj is dominated by the 

growth rate 7(k) in the regime under consideration. The 
analysis above was performed using the turbulent kinetic 
coefficients determined by Bell's turbulence that are valid if 

2 

(41) 



> K r g o , 

V V Ti J C 



where VTi - is the thermal ion velocity. A thorough discussion 
of the ef fects of a hot plasma on the short-scale modes was 
done bv lZweibel fc Everett! (|2010t ). 

We emphasize that the growth rates obtained here ac- 
count for both the CR current and the CR pressure gradient 
in the presence of short-scale Bell turbulence. In the pres- 
ence of the short-scale fluctuations, the momentum exchange 
between the CRs and the flow in the hydrodyamical regime, 
kr g o < T) , results in a ponderomotive force proportional 
to the CR current (the last terms in the mean-field momen- 
tum ecmation l28[) . As a result, there exist transverse growing 
modes with wavevectors along the initial magnetic field with 
growth rates that are proportional to the turbulent coeffi- 
cient Kt defined in Appendix [Dj In the opposite regime with 
r\ 3> 1, the fastest growth rates are for the modes having 
wavevectors nearly transverse to the initial magnetic field. 
These modes are compressible. If the short-scale Bell turbu- 
lence is absent (i.e., when the turbulent coefficients vanish), 
unstable acoustic mode s are produced, as has been studied 
earlier bv lChalovl (|l98ct ). 



6 DISCUSSION AND CONCLUSIONS 

Collisionless shocks are complex phenomena where a num- 
ber of relaxation processes are involved that redistribute the 
bulk ram kinetic energy into individual superthermal parti- 
cles and magnetic fields. Strong astrophysical shocks can 
transfer a sizable fraction of the ram energy from thermal 
particles to extremely relativistic ones and the great chal- 
lenge for modeling these shocks comes from this wide dy- 
namic range. Furthermore, if the particles have a wide range 
in momentum, the self-generated magnetic turbulence must 
have a correspondingly wide range in wavelengths. Thus far, 
no exact treatment of this process with plasma simulations 
or other methods has been possible and approximate tech- 
niques must be used. 

In this paper, we derive a mechanism for long- 
wavelength magnetic fluctuation growth in the presence of 
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a cosmic-ray current and the short-scale magnetic turbu- 
lence produced by Bell's instability. We use the anisotropic 
cosmic-ray distribution, interacting with strong Bell-like tur- 
bulence with scales below the CR gyroradii, to calculate the 
growth rate of fluctuations with wavelengths longer than the 
CR gyroradius or mean free path. As we have emphasized, 
the power in the longest wavelength turbulence is critical 
for determining the highest energy CRs a given shock can 
produce. 

The algebra needed to describe the mechanism is long 
but, schematically, this is what happens: The J cr x B/c force 
from the CR current drives the Bell short-scale instability 
at scales below the CR gyroradius. Strong, short-scale tur- 
bulence is produced if the J cr x B/c force is large enough to 
dominate the magnetic field tension in the momentum equa- 
tion @ . The short-scale turbulence influences the large-scale 
dynamics through the ponderomotive forces imposed on the 
plasma by the turbulence and the CR current. To derive the 
mean ponderomotive force one must average the momentum 
equation over the ensemble of short-scale fluctuations. When 
this is done, it is seen that there is a contribution to the 
ponderomotive force (Eqs. [D5llD8|[D9)l . In the intermediate 
regime, to first order in the amplitude of the long-wavelength 
fluctuations, this contribution is (l + Kt/-Bo)j cr x B/c, where 
B = B e z + 8b and ft = j<* ez + gf \ Tne effect of the 
short-scale turbulence appears here through the turbulent 
coefficient Kt/Bo. We find that the turbulent ponderomo- 
tive force is large enough in both the intermediate and hy- 
drodynamical regimes, and that the CR current response in 
the long-wavelength regime can no longer be neglected (see 
Section [A}- The current cannot be treated as a fixed exter- 
nal parameter, as is normally done for the short-scale Bell 
instability. 

The CR current response (Eq. IA19[) on the magnetic 
turbulence is derived here using a kinetic equation with a 
simplified collision operator. With this approximation, the 
ponderomotive force results in a long-wavelength instability 
in a way similar to Bell's instability. The angular dependence 
of the growth rate in the hydrodynamical regime depends on 
the dimensionless collision strength 77, as is shown in Fig. [21 

The method we use to average the ponderomotive and 
electromotive forces in the presence of the cosmic ray cur- 
rent is a generaliz ation of the mean fi eld method used in dy- 
namo theory (see IBrandenburg! |2009| . for a recent review). 
By introducing a parameterized relaxation time r cor for the 
triple correlations of the short-scale field fluctuations, we are 
able to obtain the mean field equations in closed form. The 
turbulent transport coefficients (Eqs. ID6I and ID7|) needed 
to determine the pondermotive force in the mean field mo- 
mentum equation are written in terms of a dimensionless 
turbulent mixing length that is dependent on r cor and the 
dimensionless turbulence amplitude Nb- 



In the intermediate regime (r] -1 < kr g o < 1), the fastest 
growing modes (i.e., Eq. I20|) appear due to the CR current 
contribution to the mean electromotive force in the aver- 
aged magnetic induction equation (|14[) . The CR contribu- 
tions are proportional to the turbulent transport coefficient 
a t . The amplification of long- wavelength magnetic fields in 
this regime is reminiscent of the large-scale magnetic field 
dyn amo model that is widely discussed in the literature (see 
e.g. iBrandenburg fc Subramanianll2005l . for a review). 

The fastest growing long-wavelength mode in the inter- 
mediate regime has non-zero helicity with a sign opposite 
to the short-wavelength Bell modes. In the hydrodynamical 
regime (kr g o < ri~ l ), the two modes with similar growth 
rates (the solid and dot-dashed curves in Fig. [T]) dominating 
the regime kr g o < 77 -1 , also have opposite helicity. While is 
is not possible to draw exact conclusions for the properties 
of the strong turbulence with our linear analysis, the growth 
rates of the long-wavelength turbulence we obtain are com- 
parable to those of Bell's instability for a reasonable range 
of parameters (compare the maxima of the solid and dashed 
curves in Fig. 0. This suggests that it may by possible to 
overcome a fundamental problem and balance the overall 
helicity with these two instabilities, but non-linear analysis 
is needed to address the issue. 

The nonlinear, long-wavelength instability requires an 
accounting of the CR response to the short-scale turbu- 
lence. Here, we have exploited the fact that the process of 
long-scale turbulence growth in the presence of short-scale 
motions and magnetic field fluctuations, has some analogy 
with t he dynamo theory reviewed recently by IBrandenburg! 
|2009l ). The key ingredient in our model, however, is the 
anisotropic CR distribution. 

Previous nonlinear work on the shor t-scal e Bel l 

instability using MH D si mulations (e.g.. iBellj 20041 ; 

Zirakash vili fc Ptuskinl 120081 ; IZirakashvili et al .1 120081 " as- 
sumed a fixed CR current as an external parameter. These 
models studied the spectral evolution in the short-scale 
range as well as the transformation of the turbulence 
through the subshock. The evolution of t he Bell modes 
downs t ream from the shock was addressed bv lPelletier et all 
l|200rJ ); lMarcowith et all l^OPd ). 

In DSA, the shock precursor containing the CR current 
has a scale length lf(p) ~ cA(p)/v 3 . Therefore, the growing 
modes are advected through the instability region on a time 
scale 

r d (p) « cA(p)/(3v 2 s ) . (42) 

To amplify a mode of wavenumber by a factor of a few, 
one needs the growth rate to satisfy 

■y(k) ■ r d > 1. (43) 

The growth rates can be estimated for a particle diffu- 
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sion model widely used in DSA (e.g., iBlandford fc Eichlerl 
where the particle mean free path is 



A(p) = r ? r 9 (p) = 3-10 12 r ? (J^J 



(44) 



Consider the energy density, E cr , for a power-law CR distri- 
bution as expected in test-particle DSA. This is the same as 
given by Eq. (|A2|) only here we normalize the CR number 
density to n cr . Then, if the power-law index is a, the CR 
energy density can be written as 



Ed = cp n CT ®(a), 



(45) 



where 



a 



$(a) = < 



a — 4 

V Po 
a-3 



pm 
PO 

for a = 4 

4 — a 



- 1 



, for a > 4 



, for a < 4, 



(46) 



and p m is the maximum CR momentum and po is the min- 
imum. It is convenient to write 



Err 6 C 



(47) 



where e cr is the fraction of the far upstream shock ram pres- 
sure transferred to the cosmic-ray energy density. Note that 
e cr varies with precursor position and with the minimum CR 
momentum, po. Then the maximum growth rate achievable 
in the model, for a mode of wavenumber k = 2ir/ A propa- 
gating at the angle cos 9 « r]^ 1 to the initial magnetic field, 
is given by 



0.1 



1.7 



1/2 



V nTy) I To 

V 2 / v. \ - 1 / 2 



-1/2 



0.01c 



(48) 



Equation (|48[) demonstrates that, for parameters typical of 
young SNRs, and for acceleration efficiencies that are high 
enough (i.e., e cr <; 0.1), long-wavelength fluctuations can 
be strongly amplified. The validity requirement given by 
Eq. (|41[) is fulfilled for most of the conditions expected in 
young galactic SNRs. 

In DSA, particle spectra far upstream in the shock pre- 
cursor close to the escape boundary can be approximated 
by Eq. (|A2|) where the minimal momentum po can be just a 
few times smaller than p m . In contrast, just upstream from 
the subshock, the minimal momentum is much lower, i.e., 
po ~ rripC. Therefore, for p m > 10 m p c, the particle mean 
free path A < 10 18 cm and fluctuations with wavelengths 
larger than A would grow with a time scale on the order of 
1,000 yr for a shock of velocity v 3 ~ 0.01c. 



Besides increasing the maximum CR proton energy a 
shock can produce, the presence of strong, long-wavelength 
magnetic field fluctuations can affec t synchrotron emission 



from relativistic e l ectrons (see, e.g ., Vink fc Laming! 120031 ; 



iBamba et al]|2005l ; lUchivama et alj|2007l ; IVinlJl200Sf ). Tern 



poral, spatial, spectral, and polarization features seen in syn- 
chrotron emission from the shells of young supernova rem- 
nants will all be modified to some extent by long-wavelength 
turbulence. The observation of these features can provide 
unique information on the properties of the long - wavelength 
fluctuations (see, e.g.. iBvkov et af]|2008l . [20091 ). In case of 
r\ > 1, and for the quasi-parallel shocks considered here, 
the growing long-wavelength magnetic field fluctuations are 
propagating obliquely to the initial magnetic field. This 
means they must be accompanied by plasma density fluc- 
tuations that are interacting with the shock front, an effect 
that potentially can be stud ied through the optica l line ob- 
servations of SNRs (see e.g.. iRavmond et alj|20101 ). 

The interpretation of optical observations might also be 
influenced, at least for the quasi-parallel shocks we studied 
here. For r] > 1, the growing, long- wavelength magnetic field 
fluctuations are propagating obliquely to the initial mag- 
netic field and they must be accompanied by plasma density 
fluctuations. These density fluctuations might produce ob- 
serv able features in optica l line observations of SNRs (see, 
e.g.. lRavmond et alj|20ich . 
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APPENDIX A: LINEAR RESPONSE OF THE 
COSMIC RAY CURRENT 

A small perturbation Sh imposed on the local mean mag- 
netic field in the form of a plane monochromatic wave results 
in a linear response of the cosmic-ray current <5j cr . General 
expressions for the linear response of a background plasma at 
re st with an arbitrary partic le distribution, were presented 
bv lKrall fc Trivelpiecej (|l973l ). Kinetic equations will be used 
here to derive the linear response of the cosmic-ray current. 
For the unperturbed CR distribution, we consider distri- 
butions tha t are typical for the diffus ive shock acceleration 
model (e.g jBlandford fc EichlerHl987f ). In this case, the mo- 
mentum distribution of accelerated particles allows a local 
approximation to a weakly inhomogeneous (on scales of the 
order of the cosmic-ray gyroradii) distribution given by 



/o(p) 



4tt 



N(p) 



1 + 



3v s p z 



cp 



(Al) 



where u s is the drift velocity of the accelerated particles 
along the mean magnetic field Bo relative to the background 
plasma ions, and N(p) is the isotropic part of the cosmic-ray 
particle momentum distribution of the power-law form 

-3 



N(p) = 



(a - 3K" 



p" 



Po < p < : 



a > 3, 



(A2) 



normalized as 



N(p)p 2 dp : 



Actually, in DSA the particle distribution is inhomoge- 
neous on the scale If ~ Ac/v 3 which is much larger than the 
particle gyroradius in the nonrelativistic shocks considered 
here. Thus, in our local analysis, the minimum momentum 
of the particle distribution po in Eq. ()A2|) is position depen- 
dent. In the far upstream region, po approaches p m making 
the distribution given by Eq. (|A2|I rather narrow. In con- 
trast, at the shock front po ~ m p c and the distribution is 
much broader. 
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The particle distribution function of energetic particles, 
f(r,p,8, ip,t), disturbed by a magnetic field, satisfies the ki- 
netic equation 

¥ +v lF + eE lF~¥ (Bo + b) - a/ = /[/] ' (A3) 

where E and b are the field amplitudes of the imposed MHD 
disturbance, and O is the momentum rotation operator, de- 
fined by 



op 



(A4) 



Al Perturbations propagating along the initial 
magnetic field 

Let 5f be the perturbation of the distribution function due 
to imposed harmonic perturbations of electromagnetic fields 
Sh and E ~ exp(ih ■ r — iuit) propagating along the initial 
magnetic field direction k = ke z = Bo/Bo- We introduce the 
CR scattering rate by magnetic fluctuations, v, assuming the 
simplest form for the CR collision operator in Eq. (|A3|) , i.e., 
/[/] = — u(5 f — Sfiso), where <5/i so is the isotropic part of the 
disturbed distribution function. Then, Sf satisfies 

- kv z )8f + - (v x Bo) ^ + vSf = 
c op 



E + 



v x Sh 



dfo 
dp ' 



(A5) 



In the case of parallel propagation, it is convenient to 
distinguish between two circular polarization modes Sh — 
5b(e x ± i&y). Then, the CR current response from Eq. (|A5j) 
takes the form 



T7E- en CI -^-v s k x 
4 uj Bq 



1 (l-t 2 )dt 
1 1 =p xt ± ia 



N(p)p 2 dp, 



(A6) 



ken v 1 eBo „ . . .., 

= — = 



eB ' 



V 



equations one obtains E = E(e x ±ie H ) = ±i — 5b(e x ±ie B ) 

kc 



We neglected E in the right hand side of Eq. (|A5f) assuming 
Vph/vs « 1. Then using j" — en CI v s from Eq. (|A1|) . an 
integration over time of Eq. (|A6[I yields 



sr = 



1'Sb / a(p)N(p)p 2 dp, 
Jo 



(A7) 



3 3 (. 1 /a\ 2 \ 3a 

CT =2^ + si i -^ + y r 1 -^*^ 



4a; 



3a 3a T 
1 *i 

2x 2 2x 3 1 



(a; + l) 2 + a 2 



(x- l) 2 + a 2 



^2(2;) = arctg ( J + arctg I 



(A8) 



For convenience, we introduce the function A(xo) de- 
fined as 



f-OO 

A(x )= / o(p)N(p)p> 
Jo 



dp 



kepo 

where xq = — — giving: 



eB 



Sj CT = g'A(x )8h 



(A9) 



(A10) 



Then, in the long-wavelength limit xo <C 1, A(xo) takes the 
form 



A(x ) 



1 



1 + i/V 



it 



In the intermediate limit , rj — > 00, Eq. 
equivalent to that of lBelll (|2004h . i.e., 



(All) 
8b has a form 



3 3 
2x 2 4x 



In 



x + 1 



4.c 



(A12) 



where 0(:r) is the Heaviside function. For the important case 
of a test particle distribution of sh ock a c celera ted CRs with 
a = 4, Eq. {A9j reduces to that of lBelll l|2004h and then 



A(xo 

_ 3 

~ 8 

37TI 



a(p)^dp 
p 



1 + 





3 
16 


Xo 


1 

T 2 


(- 


x 



In 

for xo > 1 



s 



so — 1 



(A13) 



1, for xq < 1 



In the short-wavelength limit, > 1 and A(xo) has the 
same asymptotic behavior as for finite rj as is clearly seen in 
Fig. (lAlb . 



A2 Long-wavelength CR current and pressure 

responses to oblique magnetic perturbations 

In the case of parallel propagating perturbations, the simple 
approximations for the linear response of the CR current 
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were obtained in both the intermediate and the hydrody- 
namical regimes. In the case of oblique perturbations, the 
general equations were obtained in the form of infinite se- 
quences of Bessel functions wit h cyclotron and Cherenko v 
resonant denominators (see e.g.. iKrall fc Trivelpiecelll973f ). 
Our main aim in this paper is to address the long-wavelength 
dynamics of the system with cosmic rays. If the wavelength 
exceeds the particle mean free path (the long-wavelength 
limit), or in the other words if kc/v is small, one can get 
simple analytic expressions for both the cosmic-ray current 
and the cosmic-ray pressure responses that can be used to 
obtain the wave dispersion equations. Note here that, con- 
trary to longitudinal propagating perturbations, the oblique 
perturbations induce a cosmic-ray pressure response. 
The linear response function satisfies 



dSf 
at 



+ v 



dSf 
dr 



£ 



Bo • OSf = 



-eE • ^ + ^Sh ■ Ofo - v(Sf - 5/iso). 
op £ 



(A14) 



In the long-wavelength limit, one can neglect terms of the 
order of kc/v <C 1. Now, we define a coordinate system with 
the z-axis along the unperturbed initial magnetic field, Bo, 
the polar angle 8, and the azimuthal angle <p> in the x — y 
plane. Thus, k\\ = kcosd and k± = ksind. To obtain a 
solution to Eq. (|A14[) . we neglect the electric field E on the 
right-hand-side, assuming v p h -C v s , and also assume the 
low frequency perturbation w « ct. Then, the solution to 
Eq. ()AI4|) can be presented as a superposition of spherical 
harmonics in 8 and ip with coefficients that depend on the 
particle energy, 
cc 

— 8b ■ Ofo = Si sin 8 cos tp + S 2 sin 8 sin ip, 



Si = --8b 



3v s n CI N(p) 



p " C 4"7T 

3w s n CI N(p) 



-6b 



47T 



(A15) 



Then, presenting the angular dependence of the distribution 
function as 



Sf(p) = A (p) + A 1 (p) cos 8 + 
+A2(p) sin 8 cos ip + Az(jp) sin 8 sin ip, 



(A16) 



one obtains, in the limit ckjv <C 1, the four equations for 
the coefficients: 

ikucAoM + -QAi(p) = 0, 



ik x cA (p) + -Q,A 2 {p) - OA 3 (p) = Si, 



ik y cA (p) + SlA 2 {p) + -QA 3 (p) = S 2 , 
V 



fe| t Ai(p) + k x A 2 {p) + k y A 3 (p) = 0. 



(A17) 



Using the above, the linear responses of the cosmic-ray cur- 
rent and pressure can be obtained from 



6 j = e 



j v<5/(p)d 3 ; 



5P CI = ~ / v P 8f(p)d 3 p. 



(A18) 



Substituting the solutions of Eqs. (|A17|) . into (|A16() and 
(|A18I) . and performing the integration over the particle mo- 
mentum d 3 p, we obtain finally the cosmic-ray current and 
pressure response on the imposed oblique magnetic fluctua- 
tions: 



Sf r = Y ec J o {Mp)e\\ + A 2 (p)e x + A 3 (p)e y } 



p 2 dp ■- 



fell f 1 
'— \ <^ fcii 5b + -k x 5b 



5P C ; 



4tt f° 
W 



A (p)p dp ■■ 



iB 1 



+ V^{ fc||5b+ ^ kX(5b 



(A19) 



We note that in the case of the short- wavelength limit, 
the procedure described above cannot be applied because 

— > I. In this case, one should treat cyclotron and 
Cherenkov resonances. 



APPENDIX B: CORRELATIONS IN BELL 
MODES 

To obtain the mean field dynamic equations in a closed form 
one must express the mean electromotive force £, ((vV)v), 
((bV)b), and other terms through the correlators of the 
short-scale Bell fluctuations (11* (k)6 ( g(k)) and (6* (k)6^(k)}. 
Since the maximal growth rates of the Bell modes are along 
the local mean magnetic field, it is convenient to present 
the correlations in a coordinate system with the z-axis along 
e z i = B/ I B I, where e x / ?B / are the unit vectors in the plane 
transverse to e z i . Doing this we obtain: 



(6*(k)6^(k)) = 



<b 2 (fc,')> 



xS(k x ,)S(k y >) 



\k z . 




\kA 
1 






/ 



(Bl) 
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isfies the equation 




Figure Al. The dependence of the dimcnsionlcss complex re- 
sponse function A(xq) on xq = kr g o. The imaginary parts SA(xo) 
are shown as the dotted curves, while 1 — ?RA(xo) are the solid 
curves. Left panel: the hydrodyamical case with the dimension- 
less particle mean free path t] = 2 Right panel: the limit of Bell 
(2004) corresponding to r\ — y oo. 



«(k)^(k)} 



<^(k)Mk)> = 



1 fci 
4np \k z , 



(6*(k)6 p (k)) 



(B2) 



ki (b 2 (M) 



x8{k x ,)8(k y ,) 



^p\l \k z ,\ 2 

1 \ 



\ o oo/ 



(6i(k)«p(k)) = -<«;(k)6 / j(k)>. 



(B3) 



(B4) 



APPENDIX C: MEAN-FIELD INDUCTION 
EQUATION 

To get the mean field equations we apply an averaging pro- 
cedu re that is widely used in the mean field dynamo the- 
ory (Blackman fc Field 20021 : iBrandenburg fc Subramanianl 



12005 



Brandenburg! 120091 ). The distinctive feature of our 



model is the presence of the cosmic-ray current and we shall 
derive the cosmic-ray current effect on the mean field dy- 
namics and unstable modes. The electromotive force, £, sat- 



d£ /dv , \ / db 

C m=\di xh ) + V x aF 



(CI) 



To treat d£/dt in Eq. ((CT| we obtain from Eq. (Q} the 
equation for short-scale variations of the bulk velocity, v, 
including second-order correlations: 

dv 



1 

cp- 



((j cr - en cr V) x b) + 



(C2) 



dt 

+ ^I(vxB)- (VV)v - (vV)V + 

cp 

+ -!-(V x b) x B + -!-(V x B) x b + 

Anp 4np 

+ T^(V X b) x b - J- <(V X b) X b) - 
47rp inp 

-(vV)v + <(vV)v) + — (v X b) - 
cp 

-— (v x b) + vAv, 
cp 



where j cr is the averaged cosmic- ray current. In Eq. ()C2|) . 
the short-scale density fluctuations were omitted because 
the fastest growing modes are the incompressible modes 
with wave vectors along the local mean magnetic field. The 
short-scale fluctuations of the cosmic-ray current were also 
neglected being small for kcpo/(eBo) > 1 as shown in Ap- 
pendix [X] 

The fluctuating part of the magnetic field, b, satisfies 
the equation 

9b — — 

— = V x (v x B) + V x (V x b) + 

+V x (v x b) - V x (v x b) + i/ m Ah. (C3) 

Substituting Eqs. (fC2|) and ((C3|) into Eqs. ([Cl]), ((DTj) . 
and (ID2I). one obtains: 



d£ 1 

c— = (((j cr - en cr V) x b) x b) + 

at cp N ' 



+ (v x {V x (v x B) + V x (V x b)}) - 
-<{(VV)v + (vV)V} xb) + 

+ ^ V x b ) x S ) + ((V x B) x b)} x b) + 

+E^L ((v x B) x b) + v (Av x b) + 
cp 

+v m (v x Ab) - c . 



(C4) 



The last term, £/t coi , in Eq. (|C4|I approximates the time 
relaxation of triple correlations with the tim e scale t coi 
(see, e.g.. IBrandenburg fc Subramanianl l2005h . The corre- 
lation time is typically expected to be about the turnover 
time of the turbulence. For the sake of simplicity we use 
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the same relaxation time t coi for all of the triple correla- 
tions. A more rigorous analysis that distinguishes the cor- 
relation times of different triple correlations is beyond of 
the scope of this paper and will be done separately. The 
time derivative, d£/dt, in Eq. ()C4[) suppresses the mean-field 
variations on time scales below r cor and is analogous to the 
Faraday displacem ent current in Maxwell's equations (e.g., 
lBrandenburell2009l ). This term can be omitted in the analy- 
sis of the long- wavelength modes of frequencies cjr cor <C 1 in 
the close analogy with the well-known MHD approximation. 

For Bell's instability, the terms in Eq. (|C4|) contain- 
ing the short-scale magnetic field correlators Eq. (|B1|) and 
Eq. (|B4[) (apart from a term with the cosmic-ray current) 
are ki/\k z i\ 2> 1 times smaller than the velocity correla- 
tor Eq. (|B2|) . With this approximation, the equation for the 
mean electromotive field simplifies to 



d£_ 

dt pc 

+ \ify' - en CI V y >)e y , + {fj - en CI V z >)e z , 
- (v • V x v) B z /e z , - ^Vi x B c— 

2 "7~cor 



(C5) 



where Vj_ contains transverse coordinate derivatives only. 

Using Eq. (|B2|) . the second term on the right-hand-side 
of Eq. (IC5|) can be written as 

- (v V x v) B z ,e z , = — (b 2 ) 0f, - en CI V z ,) e z , . (C6) 
pc 

Then, the mean field induction equation yields 

|=Vx(VxB) + 

+ — (b 2 ) r cor V x fif, - en CI V z A e z , + 
pc N ' v ' 



ffi - en CI V y , e 



(, 2 > 



TcorV X (V± X B) (C7) 



The long-wavelength perturbations <5j cr in Eqs. (|C5f) 
and (|D5[) were derived in a coordinate system with the axis 
e z i — B/ | B along the local mean magnetic field. To 
make a linear analysis of the dispersion relations it is conve- 
nient to transform the vector coordinates to the laboratory 
system where the z-axis is along the unperturbed magnetic 
field and the shock normal e z = B /-Bo- To first order in 
the small amplitudes of the current and field perturbations, 
the transformation yields: 

05 ~ en cr V x ,)e x , « [-g'6b x + (SfJ - en c M*)]e x , (C8) 



(j c J, - en cr V z >)e z , ss (fj - en cr V z )e z +g'Sb x e x + g'5b y e y . 

(CIO) 



APPENDIX D: THE AVERAGED MOMENTUM 
EQUATION 

To get the equation of motion averaged over the short- 
scale fluctuations in a closed form, we derive ((vV)v) and 
((bV)b) in Eq. for large-scale motions using the mean- 
field approximation as it was described in Appendix[C] This 
yields: 



i«"H = ((l')«)*(Wi). 



| <(vV)v> + — <(vV)v) 

Ob Tcor 



(bV) 



9b 



(D2) 



cp 



((((j- - en cr V) x b)V)v> 



cp 



((vV)((j cr - en cr V) X b)> 



-<({(VV)v + (vV)V}V)v)- 
-<(vV){(VV)v + (vV)V})) + 

+ 4^ ( (vV){(V x b ) x B) + ((V x B) x b)}> + 

+ T~ <({(V x b) x B + (V x B) x b}V)v) + 
Airp 

+ / ((v x B)V)v> + ^21 ((vV)(v x B)> + 
cp cp 

+i> ((AvV)v) + v ((vV)Av) , 



(D3) 



((bV)b) + — ((bV)b) = 

UL I COT 

= (({V x (v x B) + V x (V x b)}V)b) + 
+ ((bV){V x (v x B) + V x (V x b)}) + 
+u m ((AbV)b) + u m ((bV)Ab) . (D4) 



0y, ~ en cr V y >)e y , 



-g'Sby + (5j y 



C rSV y )]e y , (C9) 



Then, the equation of motion in the mean field approx- 
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imation is 



dV — — 1 

^- + (VV)V = --VP- 
dt p 

- — ((F- en cr V) x B) + -^-((V x B) x B) + 

C p i-Kp 

+ — - en cr V x i)e y , - (jf,- en CI V y >)e x ,) + 
cp 



Ct / d(fJ,-en CI V x ,) c i d(fj-en CI V y ,) ^ 
cp \ dx' dy' 



d{jf, - encrV ,/) d{fj, - en cr V z ,) 



dx' dy' 

f dV z , _ dV z , 
\ dy' x dx' 



+ 



+r cor (v • V x v) ( -^j-e x , K^~ e v' ) + 



-— ((j cr -en cr v) xb}..., (D5) 
cp 

where the turbulent transport coefficients 
f°° 1 

Kt =r cor / dfc z ,^=Vfci|fc 2 ,| (b 2 (M) , (D6) 

and 



ct=rc ° r r^'vfev^ (b2(M) (D7) 

are expressed through the correlation time r cor . 

As in Appendix [C] we transform the vector coordinates 
to the laboratory system where the z-axis is along the un- 
perturbed magnetic field and the shock normal e z = Bo/Bo- 
Then, to first-order in the small amplitudes of the current 
and field perturbations the transformation yields 

Qf, - encrVVK' « [-g'5b x + (SfJ - en CI 5V x )]e y , (D8) 
U", - en CT V v i)e x , « [-g'db y + (Sff - en CI 8V y )]e x . (D9) 
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